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00 ' Aims. We present a coherent stellar and nebular model that reproduces observations of the planetary nebula IC418. We aim to test 

whether a stellar model found to provide an optimal description of the stellar observations is able to satisfactory ionize the nebula and 
. reproduce the nebular observations, a finding that is by no mean evident. This allows us to determine all the physical parameters of 

■^1^ ' both the star and the nebula, including chemical abundances and the distance. 

Methods. We used all the observational material available (FUSE, lUE, STIS and optical spectra) to constrain the stellar atmosphere 
model performed using the CMFGEN code. The photoionization model is developed by comparing solutions provides by Cloudy _3D, 
. with results from CTIO, Lick, SPM, lUE, and ISO spectra as well as HST images. The aperture sizes and positions of the different 

observations are taken into account. More than 140 model nebular emission lines are compared to the observed intensities. The 
distance is determined using evolutionary tracks. 
, Results. We reproduce all the observations for the star and the nebula. The 3D morphology of the gas distribution is determined. 

The effective temperature of the star is 36.7±0.5 kK. Its luminosity is 7700 Lq. No clumping factor is needed to reproduce the age- 
luminosity relation. We describe an original method for determining the distance of the nebula using evolutionary tracks. The distance 
of 1.25 kpc is found to be in very good agreement with recent determination using parallax method, and the age of the nebula is 
estimated to be 1400 years. The chemical composition of both the star and the nebula are determined. Both are carbon-rich. The 
^S) ■ nebula exhibits evidence of the depletion of elements Mg, Si, S, CI (0.5 dex lower than solar), and Fe (2.9 dex lower than solar), 

^ ' which is indicative of a depletion of these elements onto grains. 

fT^ ' Conclusions. We develop the first self-consistent stellar and nebular model of a planetary nebula that reproduces all the available 

, observations ranging from IR to UV, showing that the combined approach to the modeling process leads to more restrictive constraints 

and, in principle, more trustworthy results. 
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1 . Introduction ing conditions must also be correctly simulated before compar- 

ing theoretical line ratios to observed ones. For example, if the 
Planetary nebulae (PNe) are the ultimate visible stage in the evo- observations are obtained through a slit that covers only a small 
lution of intermediate-mass stars. They consist of an extended ^^^^^-^^ ^^^^ ^ ^^^^^^^ ^^^^ theoretical 
shell photoionized by a central hot remnant (a white dwarf). ^^^^^^^^^^^ ^^^^ through such a slit. 



_i The determination of the chemical composition of PNe is im- ^, , ^ .... , , , , 

» » f • T» 11 » 1 • » » u The results from photoionization models are strongly depen- 

- - portant tor various reasons. It allows us to place important obser- , , ,^ , ..... 

. . , ... , , , »u • f • » J- » dent on the assumed spectral energy distribution (SED) of the 

vational constraints on models nucleosynthesis or intermediate- ... ^, „„„ . „ , , „ „ 

. ^/r^. . , ihnmK T» 1 11 » » ionizing Star. The SED IS often represented by a tree set or pa- 
mass Stars (e.g., Marigo_etayu003D. It also allows one to trace ° . , , ,. , • 

. „. f »u ' J- . f u- u *u * rameters in the modeling procedure. However, in some cases, 

the metallicity of the medium out of which the progenitor star , , . „ , , • , 

, -f ■ , » /c A N u u J » spectral observations of the central star exist, and allow one to 

was born, if using elements (S, Ar) whose abundances are not ^ ^ ,. , ■ , • , 

jvcj-^u » 11 1 » perform a direct analysis of the properties of the stellar atmo- 

expected to be modified in the external layers of the planetary ^ , , ^ . ^ ^ ^ . , . , 

f , /nxTN J • w c- 11 J sphere (e.g., effective temperature, gravity, chemical compost- 
nebula (PN) progenitor during its evolution. Finally, it provides - s, . , , , , 
1 u » »u J w f » 1 » J » ■ "i J u tion) by comparing the observed stellar line shapes and intensi- 
clues about the depletion of metals onto dust grams produced by ^ ° ., , „ , 
the PN progenitor or ^ t^g^ ^^^.^ those produced by a tailored stellar atmosphere model. 

Abundance determinations in planetary nebulae ai-e obtained While many detailed models of either planetary nebulae 

either by empirical methods, directly from the observed specti-a, or PNe central stars have been constructed separately, there 

or through photoionization modeling. In principle, a photoion- are only a few examples of studies combining both aspects 

ization model is expected to provide more reliable abundances, (Ercolano_etalJ |2004; Wri ght et al.| QMB and, for example, 

since these are obtained in a way that accounts for the tempera- companng the abundances determmed for the star and the neb- 
ture and ionization structure of the nebula. However, to provide 

reliable abundances, the model must reproduce the observed line In this paper, we present such a study for IC418, an el- 

ratios perfecdy, which is by no means an easy task. The observ- lipsoidal PN of high surface brightness and low excitation. Its 

rather simple morphology combined with an impressive amount 

Send offprint requests to: C. Morisset of observational material make this object an ideal test-case for 
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a detailed combined model, associating up-to-date photoioniza- 
tion and atmosphere modeling codes. 

This paper is structured as follows: a first section is devoted 
to the description of the observation material. The Section 2 de- 
scribes the stellar and nebular modeling tools, as well as the stel- 
lar and nebular modeling process and the interactions between 
the two. In the Sect. 3, we present the results of our best-fit 
model. In the Sect. 4, we discuss the expected accuracy of our 
results. In the Sect. 5, we discuss some implications of our find- 
ings, the most outstanding being related to the comparison be- 
tween stellar and nebular abundances. 



2. Observations 

2.1. Stellar observations 

The spectrum of the central star were obtained from sev- 
eral sources. The data from 900 to 1100 A were taken from 
FUSE observations obtained during the program PI 15 (PI 
J.M.Shull). A STIS/HST spectrum (dataset o52902) covers the 
range from 1190 to 1700 A and we finally combined the 
lUE spectra SWP37770, SWP37773, SWP37776, SWP37778, 
SWP37779, SWP37783, SWP37784, SWP37785, SWP37786, 
and SWP37792 to extend the UV spectrum to 1900 A. The 
existing high resolution LWP spectra are too noisy for detailed 
analysis. The photometric calibration of high resolution spectra 
were not very reliable, so we used two low resolution lUE spec- 
ti-a, SWP03178 and LWR03390, when comparing the observed 
and calculated stellar absolute flux. 

In addition to the UV spectrum, we used a high resolu- 
tion (R ~ 20000) optical spectrum of the central star obtained 
at the National Astronomical Observatory in San Pedro Martir 
(Mexico) using the REOSC echelle spectrograph. After the stan- 
dard reduction performed with ESO-MIDASQ we extracted two 
nebular spectra in windows close to the star and of same size. 
These spectra are averaged, scaled, and subtracted from the stel- 
lar spectrum. In that way, all but the strongest nebular lines were 
removed. Because of the wind velocity of ~ 500 km/s. Ha is 
broader than the nebular lines and after subtraction of the neb- 
ula, its wings were clearly seen and we use them in the following 
analysis. 

The central star of IC418 is a known irregular variable with 
two chara cteristic timescales, o ne of a few hours and one of a 
few days dHandler et al.l [19971) . The same authors report am- 
phtudes of 0.1 ma g in V and a very small amplitude in B-V. 
iMendez et al.l ( 11986) observed variations in the equivalent width 
of CIV 5808, 12A absorption lines on the shorter timescale. 
The timescale of the rapid variability is shorter than the wind- 
crossing time and is filtered by it. Patriarchi & Perinotto ( 1995|) 
analyzed several lUE spectra and concluded that the star does 
not show any detectable variability in the UV lines formed in the 
wind. These results means that the variability is principally re- 
lated to the stellar photosphere. In the analysis presented below, 
we used mainly lines formed in the wind. Based on the small 
amplitude of the changes in both brightness and color and the 
insensitivity of the P Cyg profiles, we conclude that the stel- 
lar variability does not affect significantly the model results pre- 
sented in this paper. 



' http://www.eso.org/sci/data-processing/software/esomidas/ 



2.2. Nebular observations 

Acquired in 1980-1981 with the Harvard spectral scanner, 
the CTIO 0.9m telescope observations were taken from 
Guti errez-Moreno et al. ( 1985). The Lick observa tory observa- 
tions (1991-1992) were reported by iHvung etal] (11994). High 
spectral res olution spectra exhib iting emission line profiles were 
taken from lGesicki et al.l (Il996l) . 

A deep optical echelle spectrum of IC 418 obtained in 
20 01-2002 at th e 4m Blanco CTIO telescope is described 
by ISharpee et al.l (|2004|) . who measured 807 emission lines 
of which 624 were solidly identified using EMILI software 
(ISharpee et alj|2003h . 

The HST WFPC observations were obtained from the 
database archive (proposal 8773, Hajian). Images in the Ha, 
[Nil] 6584A, and [Om] 5007A filters were used. 

IC 418 was observed in the UV by the International 
Ultraviolet (lUE) and the IR by the Infrared Space Observatory 
(ISO). We take the li ne intensities corrected for extinction from 
iPottasch et all ( 1200 4V 

The UV to optical normalization was obtained following 
iPottasch et all (|2004 . using the [O ii] 247 1/7325 A line ratio. We 
use the theoretical value given by the model (namely 0.75), tak- 
ing into account the difference between the aperture sizes of the 
optical and UV observations. The optical va lues used in this cal- 
ibration are those observed by ISharpee et al. (2004). At the end 
of the convergence process, we check that the absolute fluxes are 
reproduced. 

To connect the intensities observed by ISO to the optical in- 
tensities, we use the Balmer and Brackett hydrogen lines. These 
lines were observed with the ISO SWS spectrometer, while the 
lines of wavelength longer than 50 fim were observed by the 
LWS spectrograph. Since both spectrometers have aperture sizes 
that are larger than the entire nebula, no differential aperture cor- 
rections are needed. This is confirmed by the continuum emis- 
sion (supposed to arise from the external dust around the ionized 
nebula) being of very similar intensity when observed by both 
spectrometers in a common wavelength range (around 45 fim). 

Notice that the absolute fluxes predicted by the model for 
[On] 2471 A, the Brackett, and the Balmer lines are in a very 
good agreement with the observations, meaning that the calibra- 
tions we use to connect the IR, optical, and UV observations are 
close to the absol ute calibration . The s izes and positions of the 
apertures used by ISharpee et alj (|2004 and lHvung et alj d 19941) 
are taken into account, as de s cribed in Sect. l3.2[ 

iGutferrez-Moreno et alj ("1985') report observations in the 
Unes of [Ne v] 3425 A, [Ar v] 7006 A and He ii 4686 A. We use the 
atmosphere model to verify whether the He ii 4686A line could 
be of stellar origin. The [Ne v] 3425 A and [Arv] 7006A are not 
emitted by either the star or the nebula. 



3. The models 

3.1. Stellar model 

One of the major problems in compiling any photoionization 
model of a PN, apart from the density distribution (discussed in 
Sect. 13.2b . is the description of the ionizing spectrum. By defini- 
tion, the photons capable of ionizing the gas are never observed, 
and the simplest hypothesis involves use of a black-body distri- 
bution, defined only by an effective temperature and scaled by a 
luminosity. In an effort to understand the ionization state of the 
nebula more clearly and to resolve some discrepancies between 
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the model predictions and the observations, we use a more real- 
istic description of the ionizing flux by modelling in detail the 
atmosphere of the ionizing star. 

The stellar atmo sphere model is obtained using the 
CMFGEN code from HiUier & MiUed (11998.) . CMFGEN is a 
non-LTE line blanketed code designed for spectral analysis of 
stars with winds. The radiative transfer equation is solved in the 
comoving frame by assuming statistical and radiative equilib- 
rium. The code is capable of treating a large number of ions with 
an enormous number of individual levels. As a result, the emitted 
spectrum is reliable across a wide range of frequencies includ- 
ing the flux above 13.6 eV that is responsible for the ionization 
of the nebula. 

In contrast to other models of the star, the effective tempera- 
ture is fixed not by stellar lines but using the ionization balance 
of the nebula. Once the ionization of the nebula is reproduced, 
we then check the temperature structure of the wind using the 
lines of N III, N IV, and N V, and C III and C IV. In that way, the 
stellar model is consistent with the nebular model. It reproduces 
correctly the stellar flux bluer than the H ionization edge, which 
is not seen by the observer but strongly affects the nebula. 

The density st ructure of the win d is set by an ad-hoc velocity 
law in the form (iHilUer et alj|2003h 



y(r) = 



Star I ' / 

1 + {y phot ly core - l) exp (- (r - Rsutr) I h) 



(1) 



where Vcore represents the depth to which the model extends and 
yphoi governs the match between the hydrostatic atmosphere and 
the wind. The scale factor h is defined as 



h = 



kT{\-y) 



(2) 



where y is the mean number of electrons per atom, F = gradig is 
the ratio of the radiation pressure to the gravity acceleration, jd is 
the mean atomic weight, and is the atomic mass unit. 

The value Voo=400 km/s is set by the blue wings of all UV P 
Cyg profiles. For the C IV 1548,50A Unes, we fit the bluest part 
of the saturated absorption component (sometime called Vhiack) 
instead of the point at which the blue wing reaches the contin- 
uum called Vedge, which is usually affected by the velocity dis- 
persion in the wind and can be higher than Vbiack by 200-400 
km/s. The model with the derived value of Vk,, a lthough lower 
(by a factor of two) than the previous s tudies ( Mendez et alj 
Il988t iKudritzki et a"i]|1997t iPauldrach et al. 2004) . reproduces 
in detail a large number of lines. We attribute the discrepancy 
with other models to the higher quality of the spectra (mainly 
FUSE and HST/STIS) used in this work, which were unavail- 
able at the time of previous studies. 

When the terminal velocity has been fixed, we determine the 
value of /? = 3.0 by fitting the widfli of He II 4686A (Fig.|2]i. The 
value of h (and therefore log g) is determined by fitting the wings 
of Hy. Figure[T]shows the comparison between our model with h 
- O.O25R0, which corresponds to log g = 3.55, and the observed 
Hy wings. We use a turbulent velocity of 10 km/s during the 
modeling and then broaden the calculated spectrum with a 30 
km/s rotational velocity to fit the absorption lines observed in 
the UV spectrum. Table [T] summarize the parameters obtained. 

After fixing the wind parameters, we determine the abun- 
dances of the main elements. The mass loss rate and the He/H 
ratio is determined simultaneously from the fit to the Ha, 
He I 5876A and Hell 4686A lines (Fig. |2|. It was shown 
dPrinja etalj|2007l) that PV 1117,28A is sensitive to the wind 




Fig. 1. Comparison between observed wings of Hy (thin line) 
and the model with log g = 3.55 (thick fine) 



clumping factor. The assumed clumping factor, f = 0.1, repro- 
duces this line profile so well that we do not adjust it during the 
fitting procedure. The carbon abundance is measured from the 
fit of C III 1 178 A, C III 1247 A, and C IV 1 169 A lines (Fig|D31 
online material). The nitrogen abundance is determined from 
the NIII 1750A and NIV 1718A lines (Fig.|m] online mate- 
rial). The Oxygen abundance is set by OIV 1339,41A and OIII 
1150A, OIII 1151 A, OIII 1155A lines (Fig.|D21 online mate- 
rial). The ratios of the lines for different ionization stages of He, 
C, and N agree with the temperature determined by the nebu- 
lar response. However, the observed lines of the three ioniza- 
tion stages of oxygen imply that the temperature is ~ 2000K 
higher. On the other hand, the FUSE spectrum of IC418 shows 
strong O VI 103 1,37 A emission (Fig. lD.2l i. which is indicative of 
a super-ionizat ion effect probably affe cting all the oxygen ion- 
ization stages (iMacfarlane et al.lll994l) . The oxygen abundance 
is assumed to be a compromise between that related to the O m 
and O IV lines, because the superionization effect is understood 
to have less influence on these ions. In addition, the oxygen 
abundance affects the intensity of Hell 4686A indirectly. The 
assumed oxygen abundance is consistent with He II 4686A in- 
tensity, and we therefore assign a higher uncertainty of the oxy- 
gen abundance. 

The effect of superio nization on the wind structure is unclear 
IMacfarlane et aLld 19941) showed that the presence of X-ray emis- 
sion can alter all the ionization stages of most of the important 
ions. On the other hand, existing X-ray observations of IC418 
do not detect a point source. We know that some superionization 
exists only by the presence of O VI 1031,37A. This problem has 
no simple solution, and we leave its detailed study to a subse- 
quent paper The presence of additional ionization, which is not 
taken into account in the model, infers that all the abundances 
obtained from the model into lower limits. 

The phosphorus abundance is determined from PV 
11 17,21 A, the silicon abundance from the Si IV 1392,1402A 
and Si IV 1122A Unes, the sulfur abundance from S V 1502A 
(Fig. |D.4| i and several other FUV lines and, finally, the argon 
abundance is determined from the Ar V 1346A line (Fig. ID. 2b . 
The iron abundance is measured by fitting the iron lines between 
1250 and 1500 A. The strong lines are usually fitted reasonably 
well, but there are a large number of relatively weak lines, some 
of which are overestimated and some underestimated. We at- 
tribute these discrepancies to the incorrect oscillator strengths. 
The neon abundance cannot be determined because of the ab- 
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sence of observable lines. We define its abundance to be its solar 
value. The abundance of each element is first determined using 
the lines mentioned above and then checked by comparison with 
all other lines present in the observed spectra. The parameters of 
the model that reproduces both the nebular and the stellar spectra 
are presented in Sect. |4] 

We note that one cannot obtain a unique solution for the 
stellar luminosity by fitting only the normalized stellar spec- 
trum. The parameters describing the star are not independent. If 
the star is described as an opaque nucleus with radius Rt and 
temperature embedded in the wind then the luminosity is 
L = AncrRlT^. The effective temperature T^/f is defined to be 
the wind temperature at a radius wh ere the Rosseland o ptical 
depth, tross, reaches 2/3. It is known dSchmutz et al.llT989l) . that 
models of the same temperature T, and same transformed radius 



Rt — Rt, 



'VW25()0\^ 
, M/10-» / 



(3) 



have very similar emission-line spectra. Substituting for the ra- 
dius Rt one obtains a scaling rule for the mass-loss rate of 



Ml 
W2 



\L2 



(4) 



which implies means that an increase in the luminosity L can be 
compensated by an increase in the radius Rt and an increase in 
M. There is no easy way of removing this degeneracy except by 
knowing the distance. In Sect. 15.41 we discuss one way of de- 
riving the distance. We note that we also take into account the 
X-ray emission in the modeling of the ionizing SED. This com- 
ponent does not have an important effect on the determination 
of the main parameters of either the star or the nebula, thus this 
part is discussed in Appendix lAl 

3.2. Nebular model 

The models presented in this paper are obtained using 
Cloudy _3D (hereafter C3D), described in Morisset (200^. We 
choose to use a 3D description of the object to take into account 
the different sizes and positions of the apertures used by various 
authors, whose data we compare with the r nodel. C3D is an IDL 
(ITT) library of routines based on Cloudy dFerland et alJI 19981) . 
which produces a pseudo-3D photoionization model by combin- 
ing several radial ID models with the code Cloudy. In the present 
paper, we use version c07.02.01 of Cloudy. Each run of Cloudy 
corresponds to a given radial direction from the star, the input 
parameters being changed to reflect the geometry (e.g., inner ra- 
dius, radial position of the two shells, hydrogen density). The 3D 
nebula is then obtained by interpolation of the results of the ID 
runs, leading to a cube of physical parameters (electron temper- 
ature, ionic abundances) and emissivities of the lines of i nterest. 
More details about C3D can be found in iMorissetl (12006). 

We note that the apparent granularity of the nebula as seen 
on HST images consists only of variations of 10% in the surface 
brightness, which could be due to variations of 5% in the density. 
These small variations have no effect on the determination of the 
global properties of the nebula and we neglect them by assuming 
a smooth density distribution. 

From the images of IC418, we assume that the nebula is ax- 
isymmetric, with a 3D distribution of ionized gas that follows 
an ellipsoidal shape. The nebula should not have an inclination 
angle (that is, its axis of symmetry is in the plane of the sky). 
It is also supposed that the nebula is radiation-bounded in all 



Fig. 3. Aper t ures used to rep roduce observation s by 
ISharpee et all (|2004 (left one) and iHvung et all (1 19941) (right 
one), superimposed on the H/3 model image. 



the directions. This rather simple morphology allows us to ob- 
tain complete models with only 4 runs of Cloudy covering polar 
angles between and 90 deg. Because of the polar-axial and 
equatorial-planar symmetries of the nebula, only 1/8 of the en- 
tire object need be computed by interpolation. We use a cube 
of 40^ pixels, the entire nebula being simulated by merging the 
same cube after adequate rotations, leading to a 80"' pixel cube. 

Emission-line maps are then obtained by summing the emis- 
sivity cube in a given direction (small axis in our present case). 
To compare with the observed line intensities, the slit sizes and 
positions are also simulated. Some of the optical observations 
were obtained with apertures smaller than the entire projected 
size of the object (Hvung et al. 1994; Sharpee et al. 2004). In the 
following, we apply the same aperture sizes, orientations, and 
positions to the emission-line intensity maps when comparing 
with these author's line intensities, while the full object is taken 
into account for the UV and IR intensities. Figure [3] shows the 
apertures used. We note that from the models described below, 
we found that only 10% (3. 5%) of the H8 flux is i ntercepted by 
the optical aperture used bv ISharpee et alj (|2004|) jHvung et al.l 
1994, resp.). 

The free parameters needed to fully describe a model and 
compare it to observations are: the 3D density distribution, the 
angular size of the nebula (or its distance), the nebular abun- 
dances, the ionizing photon distribution, and the number of ion- 
izing photons emitted per unit time (Qo). The processes that de- 
termine the adopted values of these parameters are described in 
the following subsections. 

The morphology of IC418 was obtained by trial and error, 
using the HST images of H/3, [Oiii] 5007 A, and [Nii] 6584A 
and density diagnostics (11 line ratios). 

The main goal in the determination of a somewhat com- 
plex density distribution is to reproduce the central enhanced 
emission from [Oiii] 5007 A. This is reached by using a dou- 
ble ellipsoidal shell morphology, where the density is adjusted 
to reproduce the emission-line diagnostics. The density distri- 
bution cannot be recovered precisely without a high resolu- 
tion density map, as for example given by a [S 11] 4068/4076A 
map. Nevertheless, the density distribution obtained here is con- 



C. Morisset and L. Georgiev: A full stellar and nebular model for IC418 
2.5 1 1 1 1 n 2.5 1 — I 1 1 1 1 2.5 1 — I 1 1 r 



2.0 



1.5 



1.0 



0.5 



0.0 



Halpha 




6550 6560 6570 6580 



2.0 



1.5 



1.0 



0.5 - 



0.0 



Hel 5876 




5860 5870 5880 5890 
Wavelength 



2.0 




4670 4680 4690 4700 



Fig. 2. The fit to the Ha, He I 5876A, and He II 4686A lines determines the mass-loss rate and the helium abundance. Observations 
are plotted in gray, the model is in black. The absorption component blueward of He I 5876A is absent in a spectrum taken 3 days 
later The lines are saturated in the spectrum and their subtraction is incorrect. The nebular subtraction removes the weak lines 
correctly so we expect that the oversubtraction of the saturated Ha and N ii line will not affect the wings of the stellar Ha, which 
were used in the analysis. 



strained sufficiently by the overall density diagnostics and the 
surface brightness maps. 

The adopted density distribution consists of two ellipsoidal 
shells of distinct densities, sizes, and widths. Its analytical de- 
scription in a spherical coordinate system is: 

nniR < Ri„) = 0, 
riniR > Ri,„0) = no+ 

«i(0)xexp(-(:^)2)+ (5) 

n2{e) X exp(-(^^)2), 

where the angular dependencies of the parameters are as follows: 
Ri{0) = Ri xF,n-iiO),R2i0) = RixFeu-iid^nm = ni/F,„-i{0) 
and 112(0) = n2l Feii-2{ff), where 6 is the polar angle (set to be in 
the equatorial direction), and Feiiiff) is an ellipsoid magnification 
factor defined by: 

F.iUe) = ad V(sin(0))2 + {at X cos(0))2, 

where a, is an adimensional coefficient, set to be 1.3 (i=l) and 
1 .5 (i=2) for the inner and outer ellipse respectively. The numer- 
ical values of these parameters for the best-fit solution that we 
obtain are given in Table [T] 

Figure |4] shows the hydrogen density distribution in a plane 
containing the polar axis of symmetry. The full density distri- 
bution is obtained by rotation about the polar axis. We will use 
this density distribution for the model presented in the following 
sections. 

Figure ID shows H/3, [Oiii] 5007 A, and [Nii] 6584A images 
of the nebula, and the ratios of these lines, for both the HST 
observations and the model, respectively. The curves are Y=0 




-IxlO^^ -5x10^^ SxlO^'^ 1x10 



Fig. 4. Cut in the density distribution through a plane including 
the polar axis, units in cm, gray-scale in hydrogen atom/cm^. 



arcsec cuts in the surface brightness, passing through the central 
point. In the case of the HST images, the central star emission 
has been removed to allow a clearer comparison between the 
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model and the observed images. The global shape of the nebula 
is reproduced well by our model. 

The flat top of the [N ii] 6584A image is due to saturation 
effect in the observation. We artificially saturate the model image 
by taking the same ratio of the mean value at the center of the 
image (without the star) and the saturation value. The resulting 
saturated profile is shown with a thicker line, reproducing the 
observed profile. 

When the density distribution is determined, we run models 
with the stellar SED determined in Sec. 13.11 The abundances of 
the elements are then determined by fitting the most important 
lines of all the observable ions of each element. A total of 140 
emission lines are used for the model, most of them being pure 
predictions (not used for the convergence process). The com- 
plete process of convergence is achieved by fine tuning the effec- 
tive temperature until we successfully reproduce the ionization 
state of the nebula, once the H/3 absolute flux is well reproduced. 

The ISO observations show a continuum emission, which is 
commonly associated with dust. The emission reaches a max- 
imum value of 400 Jy close to 35 yum. It is beyond the scope 
of this paper to reproduce the exact shape of the dust emission. 
Nevertheless, we do include dust in the photoionization model 
to take into account its absorption when fitting the stellar contin- 
uum and determining the interstellar extinction. We check that 
the dust emission is globally of the order of what is observed. 

4. Presentation of our best-fit model 

4. 1. Parameters of the nebular and stellar models 

The parameters of the adopted model are described in Table [T] 
Table |2] lists the abundances of the adopted models (stellar and 
nebular) and values from the literature (for the nebula). These 
abundances by mass (for the model and the literature data) are 
also presented in Fig. ID. 8 1 (online only). We note that our values 
for C/H, N/H, and O/H are comparable to the highest values ob- 
tained previously, whi ch also correspond to the more recent de - 
terminations (namely IPottasch et al.ll2004t [Sterling et al.ll2007l) . 

4.2. Comparing the stellar parameters with previous studies 

Stellar models of the central star of IC418 are available i n 
the literature from'Men dez et alJ (ll988l) : lKudritzki et al.l(ll997h : 
[Pauldrach et al. (2004). The stellar parameters obtained by these 
authors and in this work are summarize in Table [3] In general, 
the parameters obtained in this paper are similar to those ob- 
tained by other authors although the values presented here have 
the advantage of being consistent with the nebula and evolution- 
ary status of the star Our value of Vco reproduces most of the 
features in the observed spectrum, especially the FUV lines that 
were not available at the time of previous modeling. In general, 
one can conclude that the parameters obtained here are compat- 
ible with the previous values and the main contribution of this 
paper is the more reliable and consistent values and the chemi- 
cal composition of the star. 

4.3. Defining tolerance and quality factor k{0) to compare 
the nebular model and the observations 

To compare the results of the model with observations, we de- 
fine errors in the observed intensities, ^, as the following: 10% 
for lines brighter than 0. 1 H/3, 20% for lines of strength between 
0.01 and 0.1 H/3, and 30% for fines lower than 0.01 H/3. To take 
into account the lowest signal-to-noise ratio of the lUE and ISO 



Table 1. Parameters of the model. 



Parameter 


Unit 


Adopted model 


TefF 


kK 


36.7 


102(0o) 


S-' 


47.5 


log(L) 


Lo 


3.88 


log(g) 




3.55 




km.s"' 


450 


Mdot 


10-*Mo.jr-' 


3.6 


Nucleus Mass 


Mo 


0.6 


Teff(X) 


kK 


290 


log(L(X)) 


Lo 


0.5 


no 


_^ 

cm 


2700 


III 


cm 


7500 


fij 


cm 


11700 


< n > 


cm"^ 


9300 


log(i?,„) 


cm 


16.09 


log(R,) 


cm 


16.16 


l0g(i?2) 


cm 


16.96 


logger I) 


cm 


i u. 


i^gV^^ 2) 




16.5 


J JC 




I 


Distance 


kpc 


1.26 


Age 


a 


1400 


mass 


Mo 


0.06 


TO 


K 


8575 






0.002 


T(0+) 


K 


8690 


T(0++) 


K 


8375 


T(N+) 


K 


8810 


T(N++) 


K 


8375 



observations and the uncertainties in the connection between the 
optical and the IR and UV domains, we add 20% to the previ- 
ously determined errors of the I R and UV lines. We also add 10% 
to the optical lines obtained by Gutierrez-Moreno et al.l (|l985h : 
Hyung et al. (1994). We note that these errors also aim to take 
into account the numerical and modeling uncertainties due to the 
limited complexity of the morphology that we adopt, the errors 
in the atomic data that we use, and the effects of the position and 
size of the apertures that we use (see the comments about rotat- 
ing apertures during the observations by Hyung et al. (1994)). 
A quality factor of the fit to an observable is defined by 

log(Omod) - log(Oob,s) 

.(O) = — , (6) 

where Omod is the value returned by the model, Oobs is the ob- 
served value, and r(0) is the accepted tolerance in dex of this 
observable. We define the tolerance to be t(0) - log(l -i- ^). 
When the absolute value of a-(O) is lower than 1, the model value 
is within the adopted tolerance of the observed value. The value 
of k{0) is also used to analyze the diagnostic line ratios. 

4.4. Line ratios and absolute fluxes 

Table |4] presents the intensity ratios obtained for our best-fit 
model and compares them with the observed values. The first 
column is the identification of the emission line. Some lines ap- 
pear more than once because they have been observed by vari- 
ous authors and through different apertures. The second column 
provides the reference of the observation, using the following 
code: for CTIO (Sharpee et al. 2004); 1 for ISO (Pottasch et al 
2004) : 2 for lUE tPottasch et al.n2004^ : 3 for Lick (.Hvung et aL 
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Fig. 5. We compare the HST images (top) and the model images (bottom) for 3 emission lines. Units are in arcsec. Gray values are 
on a linear arbitrary scale. The lines are Y=0 arcsec cuts in the surface brightness maps. 



Table 2. Abundances from the literature and for the adopted model 





He 


C 


N 


O 


Ne 


Mg 


Si 


P 


S 


CI 


Ar 


Fe 


BY NUMBER in log 


Best Stellar Model 
Best Nebular Model 


-0.60 
-0.92 


-3.30 
-3.10 


-4.24 
-4.00 


-3.42 
-3.40 


-4.00 


-4.95 


-4.29 
-4.90 


-7.11 


-5.36 
-5.35 


-7.00 


-5.29 
-5.80 


-4.74 
-7.40 



LITERATURE NEBULAR ABUNDANCES (BY NUMBER in log) 



Sterling et al. (2007) Model A 
Sterling et al. (2007) Model B 
Pottasch et al. (2004) 
Henrv et al. (2000) 
Hvuna etal. (1994) Model 
Hvung etal. (1994) ICE 
AJler et al. (19J0) 
Aller & Czvzak (1983) 
Beck etal. (1981) 
Buerger (1973) 



-0.975 
-0.850 
-1.142 
-1.066 
-1.155 
-1.155 
-0.920 
-1.030 

-0.796 



-3.22 
-3.32 
-3.21 
-3.66 
-3.52 
-3.57 
-3.22 
-3.21 

-3.70 



-3.97 
-3.97 
-4.02 
-4.07 
-4.15 
-4.08 
-4.22 
-4.13 

-4.40 



-3.29 
-3.43 
-3.46 
-3.81 
-3.74 
-3.67 
-3.37 
-3.36 

-3.60 



-3.99 
-3.75 
-4.05 
-5.02 
-5.52 
-4.22 
-4.40 
-4.13 
-4.12 
-4.50 



-5.55 
-5.16 



-6.00 



-5.40 
-5.35 
-5.36 

-5.60 
-5.56 
-5.38 
-5.38 
-5.21 
-5.00 



-7.07 
-7.11 
-6.92 

-7.00 
-7.11 
-6.85 
-7.06 

-6.82 



-5.79 
-5.81 
-5.74 

-6.26 
-5.72 
-5.34 
-5.64 
-5.74 
-5.30 



Solar (Asplund et al. 2005) 



-1.070 I -3.61 I -4.22 I -3.34 I -4.16 I -4.47 I -4.49 I -7.00 I -4.86 I -6.50 I -5.82 I -4.55 



BY MASS in log 



Best-fit stellar model 


-0.30 


-2.52 


-3.40 


-2.52 






-3.15 


-5.92 


-4.16 




-3.99 


-3.30 


Best-fit nebular model 


-0.49 


-2.20 


-3.03 


-2.37 


-2.88 


-3.75 


-3.63 




-4.02 


-5.63 


-4.37 


-5.83 


Solar (AsDlund et al. 2005) 


-0.60 


-2.66 


-3.21 


-2.27 


-2.99 


-3.22 


-3.17 


-5.28 


-3.49 


-5.09 


-4.35 


-2.93 



11994 : 4 for CTIO (iGutierrez-Moreno et al.ll985h . The third and 
fourth columns are the observed values and the results of the 
model, respectively, in unit H/3 = 100. The fifth column is k{0), 
the quality factor of the fit, as defined in Sect. 14. 31 Figure|6]shows 
all the values of k(0) versus the ionization potentials (IPs) of the 
ions producing the lines. Horizontal lines at k(0) = ±1, +2 help 
us to see the lines that are fitted within acceptable tolerance lim- 
its. 



Figure [T] shows the distribution of the values of k{0). The 
symbols are related to the emitting element. Most of the values 
of k(0) are in the [-1,1] range, that is, the value obtained by the 
model agrees with the observed one within the tolerance limits. 

We do not observe a stronger discrepancy between the model 
predictions and the observed intensities at lower intensities. We 
interpret this as an indication that the tolerances defined in 
Sect l4.3| are adequate. 
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Table 3. Stellar parameters from the literature and for the adopted model 



Reference 


log(L/Lo) 


T(kK) 


logg 


M/Mo 


d(kpc) 


M (10-^Mo/yr) 




R/Ro 


Mendez et al. (1988') 


4.2 


36. 


3.3 


0.77 


2.0 






3.25 


Kudritzki et al. (1997) 


4.3 


37 


3.3 


0.89 




2.6 


700 


3.5 


Pauldrach et al. (2004) 


4.2 


39 


3.7 


1.33 


2.0 


0.72 


800 


2.7 


This work 


3.88 


36.7 


3.55 


0.60 


1.26 


1.1° 


450 


2.17 



a) Value of smooth wind which correspond to 3.6x10 * Mo/yr with clumping factor f=0. 1. 



Table 4. Results of the Teff=36.7 kK CMFGEN model, see the text for details. 



I 


II 


III 


IV 


V 


I 


II 


III 


IV 


V 


I 


II 


III 


IV 


V 


HI 3835 





9.490 


7.841 


-1.05 


NI 8212° 





0.142 


0.002 


-16.64 


[S II ] 4070 





1.780 


2.313 


1.44 


HI 3835 


3 


6.950 


7.829 


0.45 


NIII] 1750 


2 


0.584 


0.637 


0.21 


[S II ] 4070 


3 


2.440 


3.209 


1.04 


HI 3970 





16.900 


16.511 


-0.24 


[NIII] 57.2 


1 


3.323 


2.722 


-0.59 


[S II ] 4078 





0.765 


0.750 


-0.07 


HI 3970 


3 


15.980 


16.501 


0.18 


Nil 4041 





0.012 


0.018 


1.43 


[S II ] 4078 


3 


0.880 


1.041 


0.50 


HI 4102 





24.800 


26.575 


0.73 


Nil 4176 





0.006 


0.003 


-3.43 


[S II] 6716 





2.080 


1.863 


-0.60 


HI 4102 


3 


25.240 


26.580 


0.28 


Nil 4239 





0.013 


0.012 


-0.17 


[S II] 6716 


3 


3.590 


2.501 


-1.38 


HI 4340 





44.800 


47.326 


0.58 


Nil 4435 





0.003 


0.008 


3.81 


[S II] 6716 


4 


2.090 


1.988 


-0.19 


HI 4340 


3 


47.000 


47.341 


0.04 


Nil 4530 





0.004 


0.004 


-0.53 


[S II] 6731 





4.420 


3.846 


-0.76 


HI 6563 





312.000 


291.434 


-0.71 


Nil 4552° 





0.003 


0.001 


-2.79 


[S II] 6731 


3 


6.790 


5.174 


-1.04 


HI 6563 


3 


350.000 


291.063 


-1.01 


Nil 5495° 





0.015 


0.003 


-5.99 


[S II] 6731 


4 


4.270 


4.092 


-0.16 


HI 2.62 


I 


4.718 


4.855 


0.09 


Nil 5679° 





0.067 


0.017 


-5.32 


[S III] 6312 





0.857 


0.873 


0.07 


HI 4.05 


1 


8.557 


8.679 


0.04 


NIII 4641° 





0.005 


0.000 


-38.03 


[S III] 6312 


3 


0.940 


0.925 


-0.05 


HI 7.45 


1 


2.731 


2.815 


0.09 


[OI] 5577 





0.026 


0.017 


-1.59 


[SIII] 6312 


4 


0.960 


0.865 


-0.31 


He I 4026 





2.098 


2.005 


-0.25 


[OI] 6300 





2.170 


1.137 


-3.55 


[S III] 9069 





17.800 


18.799 


0.57 


He I 4471 





4.490 


4.280 


-0.26 


[OI] 6300 


4 


2.040 


1.302 


-1.71 


[S III] 9069 


3 


14.550 


18.896 


1.43 


He I 4471 


3 


3.590 


3.288 


-0.34 


[OI] 6363 





0.759 


0.363 


-2.82 


[S III] 9532 





42.200 


46.621 


1.04 


He I 4471 


4 


4.170 


4.116 


-0.05 


[Oil] 2471 


2 


19.271 


18.374 


-0.18 


[S III] 9532 


3 


38.090 


46.861 


1.14 


He I 4713 





0.610 


0.647 


0.22 


[Oil] 3726 





123.800 


88.632 


-3.51 


[SIII] 18.6 


1 


15.204 


16.078 


0.21 


He I 5876 





13.670 


11.557 


-1.76 


[Oil] 3726 


3 


90.530 


115.323 


1.33 


[SIII] 33.4 


1 


2.617 


2.600 


-0.02 


He I 5876 


3 


10.850 


8.885 


-1.10 


[Oil] 3729 





52.340 


34.851 


-4.27 


[SIV] 10.5 


1 


1.266 


2.929 


2.49 


He I 5876 


4 


10.200 


11.102 


0.47 


[Oil] 3729 


3 


37.720 


45.228 


1 


[CI II ] 8579 





0.284 


0.183 


-1.69 


He I 6678 





3.870 


2.960 


-1.47 


[Oil] 7323 





13.740 


13.333 


-0.32 


[CI II ] 8579 


3 


0.270 


0.248 


-0.25 


He I 6678 


3 


2.550 


2.283 


-0.42 


[Oil] 7323 


3 


15.680 


18.152 


0.80 


[CIII] 9124 





0.077 


0.048 


-1.78 


He I 6678 


4 


2.750 


2.845 


0.13 


[Oil] 7332 





11.500 


10.733 


-0.72 


[CIII] 9124 


3 


0.050 


0.066 


0.80 


He I 7065 





7.040 


8.266 


0.88 


[Oil] 7332 


3 


13.380 


14.613 


0.48 


[CIII] 6162 





0.004 


0.004 


0.28 


He I 7065 


4 


5.130 


7.867 


1.63 


OI 7773 





0.075 


0.054 


-1.29 


[CI III] 5518 





0.182 


0.179 


-0.05 


He I 7281 





0.791 


0.683 


-0.56 


OI 8447° 





1.142 


0.007 


-28.20 


[CI III] 5518 


3 


0.190 


0.170 


-0.33 


[C I ] 8727 





0.033 


0.030 


-0.39 


OI 9264 





0.027 


0.020 


-1.06 


[CI III] 5518 


4 


0.450 


0.182 


-2.69 


[CII] 157. 


1 


1.100 


0.403 


-2.98 


[OIII] 4363 





0.935 


0710 


-1.05 


[CI III] 5538 





0.356 


0.356 


0.00 


CII] 2326 


2 


81.462 


88.482 


0.31 


[OIII] 4363 


3 


0.520 


0.421 


-0.63 


[CI III] 5538 


3 


0.380 


0.348 


-0.26 


CIII] 1909 


2 


27.592 


34.878 


0.89 


[OIII] 4363 


4 


0.910 


0716 


-0.71 


[CI III] 5538 


4 


1.000 


0.353 


-3.09 


CII 1335 


2 


23.212 


20.164 


-0.54 


[OIII] 4959 





72.700 


65.158 


-1.15 


[CI III] 8436 





0.006 


0.010 


1.68 


CII 1761° 


2 


2.102 


0.349 


-5.34 


[OIII] 4959 


3 


29.520 


37.300 


1.28 


[CI III] 8483 





0.011 


0.010 


-0.29 


CII 4267 





0.570 


0.452 


-0.88 


[OIII] 4959 


4 


50.100 


65.034 


1.43 


[CI III] 8504 





0.004 


0.009 


2.80 


CII 4267 


4 


0.690 


0.436 


-1.36 


[OIII] 5007 





214.000 


196.126 


-0.92 


[Arll] 6.98 


1 


5.997 


6.184 


0.09 


CII 4619 





0.011 


0.021 


2.61 


[OIII] 5007 


3 


85.870 


112.273 


1.47 


[Arlll] 5192 





0.039 


0.034 


-0.43 


CII 6580° 





0.805 


0.039 


-11.52 


[OIII] 5007 


4 


151.000 


195.754 


1.42 


[Arlll] 5192 


3 


0.030 


0.027 


-0.28 


CII 7231° 





0.169 


0.006 


-12.66 


[OIII] 51.8 


1 


15.261 


13.430 


-0.49 


[Arlll] 7135 





8.260 


7.736 


-0.36 


[NI] 5198 





0.201 


0.046 


-5.58 


[OIII] 88.3 


1 


2.407 


1.854 


-0.77 


[A- III] 7135 


3 


6.180 


5.937 


-0.15 


[NI] 5200 





0.117 


0.016 


-7.69 


on 4152 





0.018 


0.023 


088 


[A- III] 7135 


4 


6.030 


7.417 


0.79 


[Nil] 5755 





2.760 


2.599 


-0.33 


on 4341 





0.085 


0.091 


0.23 


[Arlll] 7751 





2.196 


1.867 


-0.89 


[Nil] 5755 


3 


3.910 


3.723 


-019 


on 4593° 





0.024 


0.004 


-6.37 


[A- III] 7751 


3 


0.780 


1.433 


1.81 


[Nil] 5755 


4 


2.190 


2.694 


0.79 


on 4651 





0.171 


0.172 


0.01 


[A- III] 9.00 


1 


7.334 


8.039 


0.27 


[Nil] 6548 





53.600 


53.094 


-0.10 


[Nell] 12.8 


1 


53.098 


61.004 


0.53 


[ArIV]4711 





0.003 


0.003 


0.08 


[Nil] 6548 


3 


71.230 


73.639 


0.18 


[Ne III] 3869 





3.090 


2.806 


-0.53 


[ArIV]4740 





0.004 


0.005 


1.58 


[Nil] 6548 


4 


51.300 


54.902 


0.37 


[Ne III] 3869 


3 


2.050 


1.147 


-2.21 


[Fe III] 4659 





0.027 


0.027 


-0.06 


[Nil] 6584 





162.900 


156.681 


-0.41 


[Ne III] 3869 


4 


2.630 


3.766 


1.37 


[Felll] 5271 





0.015 


0.017 


0.34 


[Nil] 6584 


3 


206.830 


217.309 


0.27 


[Ne III] 3968 





0.970 


0.846 


-0.52 


[Felll] 4881 





0.015 


0.013 


-0.52 


[Nil] 6584 


4 


151.000 


162.014 


0.39 


[Ne III] 3968 


3 


0.640 


0.346 


-1.83 


Mgll 2798 


2 


16.351 


15.697 


-016 


[Nil] 121. 


1 


0.283 


0.189 


-1 


[Nelll] 15.5 


1 


9.512 


9.276 


-0.08 


[Si II ] 34.8 


1 


0.936 


0.985 


013 



° Lines suspected to be mainly emitted by fluorescence. 



The hydrogen recombination lines are almost all reproduced 
well, with the notable exception of the ratio Ha/HyS, for which 
the observed values are 3.5 and 3.1 (depending on the observer), 
considerably higher than the predicted value of 2.9. We can sus- 
pect this to be an observational problem. 

The helium recombination lines are also mostly fitted within 
the tolerance levels. When a given line is not fitted (e.g., Hei 



5876A), this is true for only one observation of the given line, 
the other observations of the same line are reproduced well 
(/f(Hei 5876A) = -1.74,-1.09,0.47 depending on the obser- 
vation). 

The 3 most intense lines of carbon are reproduced within the 
tolerance levels. 
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The [Nn] 6548,6584A lines are reported by 3 authors, and 
the corresponding 6 observed intensities are reproduced well 
with a mean /c([Nii] 6548,6584A) of 0.13. We note the very 
small value predicted for [N i] 5200A, 8 times smaller than ob- 
served. A discussion of the electron excitation collision strengths 
of can be found in Pequignot et al. (2001), who questioned 
the values published bv lTavall (12000), but the version of Cloudy 
th at we used (na mely c07.02.01) assumes the collision strengths 
of iTavall (l2006l) , which are supposed to be 20 times greater at 
an electron temperature of about 8000 K. This line is emitted by 
the very external part of the nebula. Any small variation in the 
position of the slit or the external electron density could strongly 
affect the prediction. 

The Oi lines are globally predicted to be half of the ob- 
served values. They are also emitted close to the edge of the 
nebula, and the same considerations as for [N i] lines about the 
slit position relative to the main sites of emission should ap- 
ply. The [O ii] lines are reproduced well, with the exception of 
the [On] 3726,3729A observed bv lSharpee eTa l. (2004), which 
are observed 40% higher than the model value. On the other 
hand, the observed intensities of the same lines reported by 
iHvung etal J (119941) are 25% lower than the model values. The 
[On] 7320,7330A lines are better reproduced, with the same 
tendencies. The [Oiii] 4959,5007A lines are reported by 3 au- 
thors, the corresponding 6 intensities are reproduced with a mean 
4[Oiii] 4959, 5007 A) of 0.6. The observed value for [Oiii] 
5007A/Hy6 ranges from 0.86 to 2.1 depending on the author. 
The values predicted by the model range from 1.1 to 2.0, clearly 
showing the effect of changing the position and size of the aper- 
ture and the need to take both into account in the fitting process. 

The 5 values reported in the literature for the [Ne iii] lines 
are reproduced with a mean value for K([Ne ill]) of -0.6, while 
/f([Neii]12.8/im) ~ 0.5. 

The 10 values reported for the [S n] lines are reproduced with 
a mean k{[S ii]) of -0.12, while we found a mean a:([S hi]) of 
0.46, using 9 lines of [S iii]. The mean value of k([CI ii]) is -0.52 
(5 lines) and the mean value of k([CI iii]) is -0.24 (9 lines). For 
the Ar hnes, the mean value for A-([Ariii]) is 0.10, based on 8 
lines. 

All of these mean values for k{0) are only indicative, as some 
very faint lines can enter into the mean and dominate it. Figure|6] 
is very instructive for obtaining an idea of the global quality of 
the fit. Notice that for Mg and Si only one line of each element 
is detected. 

Some lines are permitted lines that are suspected of having 
a strong fluoresce nce contribution, and for wh ich a-(O) is much 
lower than -1 (see lEscalante & Morissetll2005L for the Nil fluo- 
rescence spectrum). They are labeled by an ° close to the wave- 
length. For example, the Niii 4641 A line is virtually not pro- 
duced by recombination of the N^^^ ion (requiring 47 eV pho- 
tons), but rather is emitted entirely by fluorescence, implying 
that the N^^ ion is responsible for the emission. 

Table ID. ll (online only) compares the observed electron den- 
sity and temperature diagnostics with the results of the model. 
The four columns correspond to I: diagnostic line ratio, II: 
observed ratio. III : model ratio, IV: quality factor k(0) (see 
Sect |4.3l ) for the model. All the density (temperature) diagnos- 
tic ratios increase with the electro n density (temperatu re resp.). 
For the optical lines, HAF refer s to lHvung etal] jl994. GCM to 
iGutferrez-Mo reno et alj (Il985h . and otherwise to Sharp ee et alj 
(2003). These results are also plotted in Fig. [8] From the 11 
density diagnostics that are available from the observed line ra- 
tios, we find very good agreement between the model and the 
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Fig. 6. Quality factors a:(0) vs. the ionization potential of the 
emitting ion. Each symbol corresponds to a given element. 



observed values. The only value of k(0) that is not in the [- 
1,1] range is [ Ar iv] 4740/47 1 1 A, perhaps indicative of a slightly 
higher density thin shell or clumps close to the star The mean 
value of A-(densitydiagnostics), without taking into account the 
[Ariv] 4740/4711 A ratio, is 0.12. 

Twenty-one diagnostic line ratios are available for the elec- 
tron temperature, and for most of them the observed ratios agree 
with the results of the model. We nevertheless note the disagree- 
ment for the low ionization species (S ii and CI ii). The extreme 
valu e of 3.1 i s obtai ned for K(T,i[On] 7320/3727 A)) for the 
Sha rpee et al.l (l2003l) obse rvation. On t he oth er hand, for the 
same diagnostic ratio using' Hvung et al.l(ll994l) observations, we 
find k(0) - -0.2. This demonstrate the difliculties of the fitting 
exercise, especially when relying on a unique set of observa- 
tions. The same apphes to the classical Te([Om] 5007/4363A) 
diag nostic, which is rath er well reproduced when comparing to 
the Sharpee et al.' ('2003') observations, but less well fitted when 
using the Hyung et al. ( 1994) or Gutierrez-Moreno et al. ( 19851) 
observations. For most of the diagnostics for which we have 
\k(0)\ > 1, there is another observation of the same diagnos- 
tic for which |a:(0)| < 1 (see Sect. 15. II for a discussion of the 
small discrepancy between the model and the observed electron 
temperature diagnostics). 

Figure ID. 51 (online only) shows that there is no systematic 
trend for k(0) as a function of either the wavelength (upper 
panel) or the critical density of the forbidden lines (lower panel). 
The first result indicates that there is no important error in the 
reddening correction applied to the emission line intensities. If 
there is a low density inter-clump medium, it could be detected 
by low critical density lines. The absence of a clear trend in the 
lower panel of Fig. |D.5| indicates that there is no detectable low 
density medium: this could mean that there are no clumps at all 
and the medium is rather homogeneous, or that any inter-clump 
medium is of a very low density (see Sects. 15.31 and 15741 for a 
discussion of the determination of the clumping factor). 

Table |5] displays the predicted absolute fluxes in the 4 lines 
used to cross-calibrate different spectral regimes, and the corre- 
sponding observed values. The agreement is extremely good. In 
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Fig. 7. Histogram of k{0) distribution. 



Table 5. Comparison of some predicted and observed absolute 
fluxes (units in 10"'° erg.cm^^.s"', except for the 6cm contin- 
uum, in Jy). The model H/3 and [O ii] 247 1 A fluxes are reddened 
using Eb-v =0 .26 and Ry - 3 . 6. The observed (reddened) values 
are taken from jPottasch et al.1 (|2004|) 



Line 


Observation 


Model 




2.7 


2.56 


[On] 2471 A 


0.24 


0.21 


Hi4.05yum 


.45 


0.59 


6 cm continuum 


1.73 


1.96 



this way we verify that our cross-calibration between the dififer- 
ent wavelength regions, described in Sect. 12.21 . is good. 

4.5. Fitting fhe lUE observations 

Figure |9] compares the observed lUE stellar spectrum with the 
stellar model used in the photoionization model. The model 
shown here is the CMFGEN star taking into account the neb- 
ular absorption and the ISM attenuation. We found that the best- 
fit extinction law to reproduce the lUE observations is obtained 




1500 2000 2500 3000 3500 

Wavelength (A) 

Fig. 9. Comparison between the lUE observations and the stellar 
SED from the CMFGEN model. 



using the iFitzpatrickl (Il999l) law with Eb-v =0.26 and Ry = 
3.6. The 2200A bump is reproduced well, as well as the EUV 
ISM attenuation. This value for Eb v is within the range ob- 
tained from the literature (from 0.14 to 0.37, see jPottasch et al.l 
[2004; .Sharpee et al.ii2004i, and references therein). 



5. Discussion 

5. 1 . Missing power 

We point out that, in contrast to the electron density, which is an 
input parameter of the model (being mainly controlled by the hy- 
drogen density), the electron temperature is an output parameter, 
as it is the result of the energy balance in the nebula, and depends 
on the ionizing spectrum and the metal abundances. Thus the 
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electron temperature diagnostics are predictions for the present 
model, because all the parameters that can have an impact on 
the temperature are fixed: the ionizing SED is constrained by the 
stellar observations and the abundances of the major coolants are 
controlled by emission lines that are already fitted. Even if the 
fits to the electron temperature diagnostic line is relatively good, 
the model presented here seems to be a little too cool. On the 
other hand, the fit to the intensities of the collisionally excited 
lines is also globally quite low (see the distribution of k(0) in 
Fig|2l). This small discrepancy between the model and observa- 
tions could be resolved by adding an extra heating process. 

We test a model with an extra heating value of lO"^^ -' x 
(nn/no) erg.cm^^^.s"' (for this we used cOS.OO version of Cloudy, 
with no - I cm"-*), to obtain electron temperature diagnostics 
closer to the observations. This leads to an increase in the mean 
electron temperature from 8575 K to 9046 K. The corresponding 
total additional power is 4 xlO^^ erg.s"', representing 30% (!) of 
the heating from the star (or the cooling of the nebula, since they 
are equal by assumption of thermal equilibrium). The collision- 
ally excited lines in this boosted model are globally closer to the 
observed values. 

5.2. Where a black-body model fails 

To test the validity of using a simple SED for the ionizing flux in- 
stead of an atmosphere model, we perform a model with a black- 
body (BB) distribution. We keep all model parameters to the val- 
ues of the model described above, except for the stellar tempera- 
ture, which is set to be 39.5 kK to produce a global ionization of 
the nebula similar to what is obtained with the CMFGEN model. 
Figures lD76l and lD77] (online only) are similar to Figs. |6] and |7]re- 
spectively, except that they involve this BB model. 

All the elements apart from Ne and Ar exhibit line intensities 
very similar to those obtained with the adopted model. For Ar 
and Ne, the ions of high IPs > 40 eV (namely Ne^^ and Ar-*^) are 
strongly overpredictedby the BB model, relative to the lower IPs 
ions from the same species. We attribute these overpredictions to 
the shape of the SED at high energy (see FigfTOl). The absorption 
features visible between 40 and 50 eV are mainly produced by 
C m, O IV, and Ne m, among other lines. 

The main difference between the BB model and the 
CMFGEN model is only visible above 40eV. For energies below 
this limit, the two SEDs are very similar. Thus, in the absence 
of any atmosphere model, the Planck hypothesis is a good ap- 
proximation to the ionizing flux between 10 and 40 eV, but there 
is a significant shift in the effective temperature (here a Planck 
function at 39.5 kK is similar to an atmosphere model close to 
36.7 kK). 

We note that i Hvung et al.l (1 1 9941) obtained very low abun- 
dances for Ne and Ar using a photoionization model of IC418. 
We suspect that the SED that they used is responsible for this 
discrepancy with our results. 

5.3. Clumping factor and distance determination 

None of the geometrical size of the nebula, the distance, nor 
the absolute luminosity of the central star are well-constrained. 
Multiplying the stellar luminosity by a factor of Ki but instead 
all the sizes of the nebula (inner radius, positions and sizes of 

1 /3 

the 2 shells) hy will lead to the same position of the recom- 
bination front and the same appearance (because the Stromgren 

1/3 

radius is proportional to Qq ). If the distance to the observer is 
also multiplied by K^^^, the angular size remains constant, while 



the absolute luminosity (Hy8 for example) is multiplied by the 
same factor . On the other hand, changing Qo and the sizes 

1/3 

by factors of Ki and , respectively, leads to increase in the 

ionization parameter of U - Qo/(4.k.R^.Nh.c) oc K^^^. This 
will lead to an increase in the intensities of the lines emitted by 
high IP ions (e.g., Oni, Ariv, Nem) and a decrease in the in- 
tensities of the low IP lines (e.g., S n). This behaviour is very 
similar to the reaction of the nebula when the stellar temperature 
increases. We use the H/3 absolute flux to resolve the degeneracy 
and to fine-tune the effective temperature of the star (to within an 
uncertainty level of 500 K, the first order value being obtained 
from the stellar spectrum fit) by reproducing the ionization state 
of the nebula. 

Another factor that can strongly increase the uncertainty in 
the morphology (and distance) determination is the filling factor. 
This factor can be divided into 2 parts: a morphological filling 
factor //m and a clumping factor ffc. The first one represents 
the ratio of the volume of the gaseous nebula to the volume of the 
total extension of the nebula. In the simple case of a thin shell of 
external radius R„eh surrounding an empty cavity of radius Rem. , 
we have //m = l-(-§^)^. The use of a detailed description of the 
nebular density structure fixes this morphological filling factor 
On the other hand, the clumping filling factor aims to take into 
account any possible small-scale structure of the nebula, where 
the gas would be concentrated into small clumps. The physical 
process that could lead to this situation is unclear; the real situ- 
ation is certainly not clumps in the middle of a vacuum, but this 
filling factor is often used as the first moment of a more com- 
plex distribution. This is the goal of the "filling factor" keyword 
used in Cloudy. The global action of this factor is to artificially 
increase the geometrical size of the nebula. There is virtually no 
difference between observables of a model with Qq, S and ffc 
and a model with x Qo, K2 x S and f fdKi, where Qo is 
the number of ionizing photons emitted by the central source, 
5 is a parameter describing all the sizes of the nebula {R\ and 
R2 for example in Eq. |5]l, and K2 is a coefficient greater than 1 . 
All the line intensities (absolute and relative), the angular sizes, 
and the images are exactly the same in both models. The nebula 
would be at a distance D in the first case and K2y. D m the sec- 
ond. Once the existence of this clumping factor is accepted, then 
there is no way of determining the distance to the nebula, nor the 
absolute luminosity o f th e star, in contrast to t he claims of either 
iMonteiro et al] (l2005l) or lGesicki etall (l2006h . 

Using a filling factor of 1.0 as a first guess (i.e. no clumps), 
we determine the following parameters: the Hy6 flux value (with- 
out reddening) is predicted to be 7.2 x 10"'° erg cm"^ s"'. The 
distance is estimated to be K2 x 1.26 kpc. The ionized hydro- 
gen mass IS Kl X 0.057 solar masses. With an equatorial radius 
of X 1.1 X lO'^cm and considering an expansion velocity of 
30 km/s (see Sect. [B), the age of the nebula is K2 x 1400 years. 
The stellar luminosity is x 7685 solar luminosities and the 
corresponding number of ionizing photons is x lO^'^'^s 

Because of the presence of stellar wind, it is impossible to 
obtain the absolute luminosity of the star and so likewise impos- 
sible to determine its spectroscopic distance. 

Th e distance to IC418 was determined by iGuzman et al.l 
(I2009h . using parallax expansion between two VLA observa- 
tions with more than 20 years of time delay. They determined 
a distance of 1.3 ±0.4kpc. This value leads to a determination 
of the parameter K2 close to unity and the other -dependant 
parameters as given in Table [T] 
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Fig. 10. This is the comparison between the CMFGEN model and a black body at 39.5 kK. 



In the following section, we present an independent (and 
more indirect) way of determining the distance to the PN. 



5.4. Using evolutionary tracks to determine tine distance 

Once the effecti ve temperature is known and us ing the evolution - 
ary tracks from'Vassili adis & Wood ( 11994 or lBloeckeij (1 19951) . 
a simple relation between the stellar luminosity and the evo- 
lutionary age of the nebula can be obtained. For Tgff close to 
36.7 kK, the relation is l og(age)= 17.34 - 3.73 x log(L/L0) 
from lVassiliadis & WoodI (fl994) and log(age)= 20.6 - 4.53 x 
log(L/Lo) from 'BloeckeiJ (11995). 

On the other hand, from Sec l5.3l we determined the age to be 



K2 X 1400 years and the luminosity to be K2 x 7685 solar lumi- 
nosities. We then have: age(years)= 1400 x yfL/TEES, similar to 
log(age)= 1.20 + 0.5 x log(L/L0). 

The intersection of this relation with those determined from 
evolutionary tracks lead to a determination of the age, the lu- 
minosity, and all the parameters dependi ng on K^. A value 
of K2 - 0.92 (0.97) i s found when using IVassihadis & WoodI 
( 1 19941) dBloeckej 19951 respectively), leading to a clumping fac- 
tor very close to 1 . 

The results presented in this section strongly depend on the 
evolutionary tracks on the one hand, but on the other hand also 
on the determination of the dynamical age of the nebula. This lat- 
est is supposed to be the ratio of the size of the nebula to the ex- 
pansion velocity. This implicitly assumes that the expansion ve- 
locity is constant during the expansio n of the nebula , which is the 
simplest hypothesis. Observations bv lRicher et al.l (I2OO8 ) seems 
to indicate that the expansion is accelerating, confirming theoret- 
ical predi ctions on the RAM pressure acceleration of the AGE 
envelope (iGarcia-Segura et al. 2006, and references therein). In 
this case, the age that we determine from the line profiles is only 
an upper limit to the true dynamical age of the nebula, and this 
will induce a shift in the position of IC418 in the evolutionary 
tracks. For example, if we consider that the expansion velocity 
is only half the value that we used in the previous section, to 
take into account the acceleration, the value obtained for K2 is 
reduced by 0.03 dex. We note that this decrease in the mean ex- 



pansion velocity leads to a value of K2 lower than 1, and the 
corresponding filling factor will he I/K2 > 1. 

The method presented in this section for determining a dis- 
tance from evolutionary tracks clearly depends on the reliabil- 
ity of these latest. The results a r e com par able when the t racks 
from either IVassihadis & WoodI (11994 or lBloeckeil ([1995 ) are 
used. The derived filling factor is close to unity, and the dis- 
tan ce is comparable to t he one obtained by the parallax method 
by iGuzman et alj (l2009i) . which all implies that the evolution- 
ary tracks are trustable, at least in the part of the HR-diagram 
occupied by IC418. 

5.5. Comparing the stellar and the nebular abundances 

The nebula is found to be both carbon- and nitrogen-rich (0.5 
and 0.2 dex above the solar value, respectively) while the oxygen 
abundance is close to solar (0.05 dex lower). The neon and argon 
abundances are close to the solar values. Magnesium, silicon, 
sulfur and chlorine are underabundant by 0.5 dex relative to their 
solar values, while iron is deficient by 2.9 dex, suggesting the 
depletion of those elements in dust grains (whose presence is 
supported by the IR continuum emission). The nebular helium 
abundance is 0.2 dex greater than its solar value. 

The star is also found to be carbon-rich. The values obtained 
for the ratio C/O (by mass) are comparable: 1 . and 1 .5 for the 
star and the nebula respectively. The uncertainty in C/H is quite 
large in the case of the nebular value. If we only consider the 
recombination lines of C, we would obtain higher value for C/H, 
while the forbidden lines would imply a lower value. A closer 
agreement between the stellar and nebular values of C/O is then 
achieved when using only the forbidden lines. 

We developed a test model by adopting the same abundance 
for the star and the nebula (equal to the mean value of the two 
abundances), and the model was unable to fit the observations. 
This implies that the differences between the star and nebula are 
greater than the tolerance level we could adopt in the fit of the 
observations. 

We also highlight that the presence of the X-ray component 
(SectjAj is not taken into account in the stellar model process. 
The true ionization state of the irradiated atmosphere can be 



C. Morisset and L. Georgiev: A full stellar and nebular model for IC418 



13 



higher than the one predicted by the model. The effect on the 
abundances determined for the star is hard to estimate, but would 
be in the sense of underestimating of the true abundances. 

The nebular composition should reflect that of the star dur- 
ing the epoch it was ejected. Therefore, the observed differences 
in the compositions should reflect the evolution of the star af- 
ter the ejection. First of all, the He abundance in the star is 
higher than that of the nebula. One possible reason for this is 
the stellar wind. The star lose part of its hydrogen envelope in- 
creasing its relative abundance of helium. Carbon is enriched in 
both the star and the nebula, which means that the nebula was 
ejected after the third dredge-up (TDU) when the carbon-rich 
material is exposed at the surface. The nitrogen abundance is 
slightly increased but N/O ~ 0. 1 . This low value, together with 
the normal oxygen composition means that no hot bottom burn- 
ing (HBB) occurred (Karakas et al. 2009). The small increase of 
N can be attr ibuted to TDU and mixing from a incomplete H 
burning shell jCristallo et al.ll2009.) . Thus, CNO abundances are 
in accordance with the progenitor mass being less than 4.0 (the 
limit for HBB). In the previous section, we calculated the age 
of the nebula and the luminosity of the star from evolutionary 
tracks. The solution obtained implies proge nitor masses of 1 .75 
Mq an d 3.0 Mp calculated w ith the tracks of lVassiliadis & WoodI 
(1 19941) and lBloeckeij (1 19951) . respectively. Thus, these progenitor 
masses are in agreement with those predicted for models that re- 
produce the observed CNO abundances. 

6. Conclusion 

We have presented a combined stellar and nebular model of the 
planetary nebula IC4 1 8 that reproduces all the observational ma- 
terial available: optical and UV stellar observations (lines and 
continuum) and IR, and deep optical and UV emission-line ob- 
servations as well as HST images for the nebula. A total of 140 
emission lines are considered, most of them being predicted en- 
tirely by the model (i.e., are not used for the convergence pro- 
cess). 

The use of a 3D photoionization model for the nebula is not 
essential for describing the morphology, since IC418 is close 
to spherical, but it is of primary importance when accounting 
for the effect of the position, size, and orientation of the aper- 
tures used in the different observations. The apertures are used 
as masks over the emission-line surface brightness maps to com- 
pute the line intensities. The apparent discrepancies between the 
intensities of a given line reported by different authors are some- 
times reproduced well when these slit effects are taken into ac- 
count. 

The agreement between the stellar atmosphere model per- 
formed using the CMFGEN code and the photoionization mod- 
els obtained by Cloudy _3D is impressive: the high ionization po- 
tential ions (Ne^^ and Ar^+), whose lines are overpredicted when 
using a simple Planck function as the SED, are perfectly repro- 
duced when a stellar atmosphere model is used. This reflects the 
effect of line blanketing in the > 40 e V region, which drastically 
reduces the ionic fraction of these ions. 

The possible presence of clumps in the nebula leads to a de- 
generacy between the distance, the luminosity, and the clumping 
factor It is a priori impossible to resolve this discrepancy in the 
case of IC418 because no spectroscopic distance can be obtained 
because of another degeneracy between the stellar luminosity 
and the mass-loss rate. Nevertheless, using the distance deter- 
mined by the parallax method by Gu zman et a l. (2009), we can 
resolve the degeneracy and conclude that there are no clumps in 
the nebula. 



We found that evolutionary tracks allows us to complete 
another more indirect determination of the distance, close to 
1.25 kpc and the age of 1400 years. This distance is in very good 
agreement with the one obtained using the parallax method. 

High spectral resolution nebular emission-line profiles are 
reproduced by a velocity expansion law that strongly increases 
with the distance from the center of the nebula. 

Reliable abundances are determined for the nebula and the 
star. We found the nebula to be carbon- and nitrogen-rich while 
the oxygen abundance is close to solar, as are those of neon 
and argon. Magnesium, silicon and iron are highly underabun- 
dant relative to their solar values, which implies that there has 
been some depletion of those elements into dust grains. Sulfur 
and chlorine are also deficient. The nebular helium abundance is 
0.2 dex higher than the solar value. The star is also found to be 
carbon-rich. 
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Appendix A: Presence of an X-ray component 

The UV spectra exhibit a strong stellar OVI 1031,1037A line (Fig. (02) that 
cannot be reproduced by a ~ 40 kK stellar atmosphere model. This line is classi- 
cally explained by the Auger ionization due to X-ray emission produced by wind 
instabilities. A weak X-ray emission is confirmed by Chandra observation^. 
The Chandra data reduction has been performed using the Ciao V3.4 package 
CALDB (updated in May 2007) following standai'd procedures. Images in the 
0.5-1.5 keV, 1.5-2.5 keV, and 2.5-8 keV bands were generated using "dmcopy" 
and "csmooth" Ciao procedures. Counts are only significantly detected in the 
softer energy band image (0.5-1.5 keV), indicating that the source is intrinsically 
soft. 

The limited count detection (less than 40 photons) does not allow a robust 
spectral analysis. The soft energy flux (i.e. 0.5-2 keV band) is estimated by as- 
suming an emission model of a black-body with a temperature of ~ 0.03 keV and 
taking into account a neutral absorption of Nh = 2 x 10-' cm"-. The unabsorbed 
measured flux of the source in the 0.5-2 keV band is 1 X 10"''*erg/s/cm2. 

To take this X-ray emission into account, we add a black-body emission at 
300 kK to the stellar flux. We do not have any constraints on the shape of the 
X-ray SED emission, so we used the Planck function as the simplest SED. 

This soft X-ray component does not affect the determination of the nebular 
parameters, such as the abundances. It does have a small influence on the high 
ionization potential (IP) emission lines ([Nelll] and [ArlV] ), which are sUghtly 
underestimated by the model if the X-ray component is not taken into account. 
We adjust the luminosity of the soft X-ray component, so that these high IP 
emission hnes are reproduced. The required bolometric luminosity is 3 L© (using 
1.25 kpc for the distance), which corresponds to a flux in the 0.5-2 keV band of 
4 X 10""erg/s/cm^, much lower than the observed X-ray flux. A hotter X-ray 
source would have a lower luminosity to conserve the soft X-ray flux and no 
influence on these lines. In contrast, a cooler X-ray source would require a higher 
luminosity, leading to intensities for these high IP lines that are higher than the 
observed values. 



Appendix B: Line profiles 

IGesicki et alj 11996) published observations of hi gh spectral resolution of 1C4 1 8 
for the HjS, [0 ni] 5007 A, and [N ii] 6584A Unes. lSharpee et alj <2004[) also pre- 
sented high resolution observations for various ions. In both cases, the H i pro- 
files are broad, with half width at half maximum (HWHM) close to 18km/s, the 
[O III] lines are narrower, with a HWHM close to 9 km/s, while the low density 
species, such as [Nii], show a double-peaked profile, with an inter-peak spacing 
of 20 km/s and each peak having a HWHM comparable to the one of the single 
[0 III] line. 

We compute th e line profiles and compare them to the observations of 
IGesicki et aIHl99 a) using the same aperture size and position they used, namely 
a centered square aperture of 3 arcsec size. Despite the density law being deter- 
mined by these authors in no way reproduces the observed surface brightness, we 
can fit an expansion velocity law similar to the one that they adopt. The resulting 
profiles reproduce the main properties of the observed profiles (see Fig. IB.l) . 
The velocity law we adopt here isVccR^ with a maximum of 30 km/s at the 
outer edge of the nebula. It is beyond the scope of this paper to make a complete 
dynamical model of 1C418. We only check that the morphology we used is not 
trivially ruled out by the fine profile observations. The values for the velocities 
of O** and that we found are compatible with the radiation-hydrodynamic 
models presented bv lSchonbemer et alj 4200^ . 

Appendix C: Chlorine collision strengths and 
electron temperature diagnostic 

The c ollision strengths for Cl^^ used in Cloudy are taken from 'Wils on & Belj 
12002 ). Note that Cloudy takes into account the correction by a factor of 10 
for the collision strength for 'Dt-'Sq as discussed in Keenan et al. (2003). We 
test these values against the previous values from Krueger & Czyzak , ( 1970j) 
and found that these latest are in closer agreement with the observations. On 
one side, the values for [Clii] and [Cliii] are fitted when using values from 
[Krueger & Czyzak , while it is impossible to fit both ions w ith the same CI abun- 
dance when using collision strengths from lWilson & Belli . In this case, the [CI ii] 
lines are predicted to be too high, while the [CI iii] lines are fitted well. 

On the other hand, the elecfi'on temperature diagnostic [CI n] 6 1 62/9 1 24A is 
found to be predicted co herent with the [S n ] and one of the two [O n] diagnostic 
when using the data of iKrueger & Czvzakl . For these low ionization lines, the 
model overpredicts the ratios with ^([Sn], [Clll]) ~ 1.5. Using the values from 




-40 -2C a 2D 40 



Fig . B.l. Emission- line profiles: models (lines) and observations 
CGesicki et al.lll996i, symbols), for Hy6 (black solid line and dia- 
monds), [O m] 5007A (blue dot-dashed line and triangles), and 
[Nil] 65 84 A (red dashed line and squares). Intensities are scaled 
so that the maximum of each profile reaches 1.0. 



Table D.l. Diagnostic line ratios, see text for the columns defi- 
nitions. 



I 


II 


III 


IV 


R0 2 2471/7325-H 


0.7635 


0.7635 


-0.00 


NeS 2 6731/6716 


2.1250 


2.0641 


-0.12 


NeS 2 6731/6716 HAF 


I.89I4 


2.0693 


0.25 


NeS 2 6731/6716 GCM 


2.0431 


2.0588 


0.02 


Ne 2 3726/3729 


2.3653 


2.5433 


0.55 


Ne 2 3726/3729 HAF 


2.4000 


2.5500 


0.24 


NeS 3 18.7/33.5 


5.8102 


6.1838 


0.15 


Ne CI 3 5538/5518 


1.9560 


1.9852 


0.04 


Ne CI 3 5538/5518 HAF 


2.0000 


2.0465 


0.05 


Ne CI 3 5538/5518 GCIVI 


2.2222 


1.9420 


-0.30 


Ne 3 51.9/88.4 


6.3413 


7.2426 


0.33 


NeAr 4 4740/47 11 


1.2000 


1.7769 


1.11 


TeS 2 4068-t-/6731-i- 


0.3915 


0.5367 


1.56 


TeS 2 4068+/6731-l-HAF 


0.3198 


0.5540 


1.92 


TeC12 6162/9124-t- 


0.0108 


0.0182 


1.60 


Te O 2 7320+/3727+ 


0.1433 


0.1949 


3.10 


Te 2 7320-H/3727-(-HAF 


0.2266 


0.2042 


-0.55 


Te N 2 5755/6584-1- 


0.0127 


0.0124 


-0.15 


Te N 2 5755/6584-h HAF 


0.0141 


0.0128 


-0.32 


TeS 3 6312/9531-h 


0.0143 


0.0133 


-0.25 


TeS 3 6312/9531-hHAF 


0.0179 


0.0141 


-0.67 


Te S 3 6312/18.7 


0.0564 


0.0543 


-0.11 


TeAr 3 7136/8.99 


1.1262 


0.9620 


-0.43 


Te Ar3 5192/7136H- 


0.0037 


0.0036 


-0.10 


Te Ar3 5192/7136-H HAF 


0.0043 


0.0037 


-0.39 


Te N 3 1750/57.4 


0.1757 


0.2336 


0.58 


Te 3 4363/5007 -h 


0.0033 


0.0027 


-0.68 


Te O 3 4363/5007 -H HAF 


0.0045 


0.0028 


-1.32 


Te O 3 4363/5007 -H GCIVI 


0.0045 


0.0027 


-1.40 


Te 3 5007-H/88.3 


119.1307 


140.8424 


0.49 


Te 3 5007-^/88.3 HAF 


47.9473 


80.5872 


1.45 


TeNe 3 3869-^/15.5 


0.4268 


0.3934 


-0.23 


TeNe 3 3869-^/15.5 HAF 


0.2828 


0.1607 


-1.47 



^ ID: 7440, exposure time of 50 ks, obtained the 12th of December 
2006 



IWilson & Bell leads to a prediction of the [CI n] 6162/9124A ratio that is lower 
than what is observed, with k([C\ ii]) 1.2. 
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Fig.D.l. The fits to the nitrogen and FelV hnes. Note the nebular emission line CII 1761 A. 

Appendix D: Online figures and tables 
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Fig.D.3. The fits to the carbon lines. NIV 1183,85A and N V 1188A hnes are also marked. 
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Fig. D.4. Phosphorus and silicon hnes. Note the strong absorption hnes. 
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Fig. D.5. Quality factors k{0) versus the wavelength (upper panel) and critical density (bottom panel). 
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Fig.D.6. Same as Fig.|6] but for a black-body model. 
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Fig.D.7. Same as Fig.|7] but for black-body model. 
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Fig. D.8. Nebular abundances by mass from the literature (small crosses, same as in Table |2]l and from the adopted model (large 
stars represent the stell ar abundance, and lar ge diamond the nebular abundances). Units are in log(abundance by number) relative 
to the solar values from lAsplund et all (l2005h . 



